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We analyze theoretically the measurement of the mean output work and its fluctuations in a recently proposed 
optomechanical quantum heat engine [K. Zhang et al Phys. Rev. Lett. 112, 150602 (2014)]. After showing 
that this work can be evaluated by a continuous measurements of the intracavity photon number we discuss 
both dispersive and absorptive measurement schemes and analyze their back-action effects on the efficiency of 
the engine. Both measurements are found to reduce the efficiency of the engine, but their back-action is both 
qualitatively and quantitatively different. For dispersive measurements the efficiency decreases as a result of 
the mixing of photonic and phononic excitations, while for absorptive measurements, its reduction results from 
photon losses due to the interaction with the quantum probe. 
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I. INTRODUCTION 

The thermodynamic description of quantum heat engines 
(QHE) has been discussed at least since the early days of laser 
physics ill] and has recently attracted much interest lH-H], 
in part because the increased control achievable over micro¬ 
scopic and mesoscopic systems opens promisi^ new avenues 
of theoretical and experimental investigation llbHlll] . 

QHE can exhibit intriguing properties, including their po¬ 
tential to outperform their classical analogues. For example, 
it has been shown that a quantum photo-Carnot engine can 
extract work from a single reservoir if the latter has built- 
in quantum coherence ifl^ . and its power can be increased 
by noise-induced coherence Q. In a different situation, a 
trapped ion based quantum engine operating on an Otto cycle 
was shown theoretically to break the Carnot efficiency limit in 
the presence of a squeezed reservoir IB. 

The definition of thermodynamical quantities in the quan- 
tum context presents however conceptual challenges |[14 | - [f^ . 
and much attention has been devoted to the proper defini¬ 
tion and the quantum statistical properties of quantities such 
as heat, work and entropy In closed quantum sys¬ 

tems work may be defined in terms of a two-time measure¬ 
ment scheme ll29l - [^ or, in a recently proposed alternative ap¬ 
proach, of a single projective measurement if^ . However the 
situation is less clear for open quantum systems, where there 
are still open questions regarding the definition and and exper¬ 
imental measurements of work and heat due to 

the lack of energy conservation in the reservoir(s). In this con¬ 
text quantum stochastic thermodynamics like its classi¬ 
cal counterpart [s^, offers an interesting framework to 
discuss thermodynamic properties and simulate numerically 
the system behavior. 

Optomechanical systems are prime candidates to investi¬ 
gate the properties of QHE. Thanks in particular to remarkable 
advances in nanofabrication they have witnessed rapid devel¬ 
opments in the last decade and can now operate routinely deep 
in the quantum regime, with broad potential for applications 
in quantum technology ll^ . Recently, three of us proposed 


and analyzed theoretically an optomechanical QHE based on 
an Otto cycle j^ . In this system the intracavity field of an 
optomechanical resonator interacts coherently with a single 
mode of vibration of a mechanical resonator. The basic idea 
behind the thermodynamic cycle is that depending on the de¬ 
tuning between the driving optical field and the resonator the 
nature of the normal modes of the system (polaritons) can 
be changed from photon-like to phonon-like, with associated 
coupling to thermal reservoirs of different temperatures (cold 
for the photons and warmer for the phonons.) The exquisite 
experimental control that can be achieved in optomechanics 
suggests that such a QHE may be a good candidate to imple¬ 
ment a measurement of the work output. 

We consider that specific system to discuss several aspects 
of the work that can be extracted from QHE. Particular em¬ 
phasis is placed on the quantum measurement of the work 
and its fiuctuations, and also on the back-action of its mea¬ 
surement on the efficiency of the system. We first show that 
for that specific QHE the work can be evaluated from mea¬ 
surements of the intracavity photon number. We consider and 
contrast an absorptive and a dispersive measurement scheme, 
both involving passing a stream of two-state atoms through the 
resonator. The former situation results in projective measure¬ 
ments of the photon number, and the associated coupling be¬ 
tween the normal modes of the optomechanical system, while 
the latter corresponds to the addition of an additional energy 
dissipation channel for the photons. We numerically deter¬ 
mine the mean work and its variance over the entire thermo¬ 
dynamical cycle for both measurement schemes and use these 
results to evaluate the measurement back-action in the thermo¬ 
dynamic cycle. Our analysis is carried out within the frame¬ 
work of quantum stochastic thermodynamics, with the mea¬ 
sured work evaluated via continuous detection of the mean 
photon number in the cavity which is responsible for the radi¬ 
ation pressure acting on the mechanical resonator. 

The paper is organized as follows. Section II briefiy reviews 
the optomechanical QHE of Ref. ll^ and the main features of 
the Otto cycle. In particular we draw attention to the fact that 
the two normal modes of the system undergo two distinct ther- 
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modynamic cycles. This will be important to keep in mind in 
the context of the measurement back-action of dispersive con¬ 
tinuous measurements. Section III defines the work output of 
the engine, using the conceptually simpler case of classical 
measurements to show the relationship between the extracted 
work and the mean intracavity photon number. This result is 
used to justify the use of continuous measurements of the in¬ 
tracavity photon number operator to determine the expectation 
value of the work and its fiuctuations in the quantum regime. 
Section IV introduces two specific types of continuous mea¬ 
surements that involve either the dispersive or the absorptive 
interaction between the cavity mode and a stream of two-state 
systems. Information on the intracavity field is then inferred 
from measurements of the state of the atoms as they exit the 
optomechanical resonator. Section V summarizes the results 
of numerical simulations obtained by a standard quantum tra¬ 
jectory approach to the solution of the stochastic master equa¬ 
tions describing the continuous measurements. Finally Sec¬ 
tion VI is a summary and outlook. 

II. OPTOMECHANICAL OTTO CYCLE 

This section briefly reviews the main features of the op¬ 
tomechanical quantum heat engine of Ref. JUl]. We consider 
a standard optomechanical setup with a cavity mode of fre¬ 
quency ojc and damping rate k coupled via radiation pressure 
to a single oscillation mode of a mechanical resonator of fre¬ 
quency (jJm and damping rate y. The cavity is driven by an op¬ 
tical pump field of strength Ofin and frequency ojp. We assume 
that the system reaches a steady state with mean intracavity 
field amplitude a ^ ain/A, where 

\^\ = \(Up-(Uc-2Pg\:^ K ( 1 ) 

is the detuning between the pump and cavity fields corrected 
for the equilibrium position of the mechanical oscillator, p = 
-ga^/ojm, and g is the single-photon optomechanical coupling 
constant. 

Denoting the small fiuctuations of the photon and phonon 
modes around the steady state by the bosonic annihilation op¬ 
erators a and b respectively, the Hamiltonian for these fiuctu¬ 
ations is 

Hab = -bAa^a fiojmb^b fiG{a^ + a){b^ + b), (2) 

where we have introduced the linearized optomechanical cou¬ 
pling strength G = ag. 

We focus on the red detuned regime A < 0, which in gen¬ 
eral leads to stable dynamics for small damping and per¬ 
form a Bogoliubov transformation to diagonalize the Hamil¬ 
tonian ([2]) in terms of normal modes (polaritons) described by 
the bosonic annihilation operators A and B. Ignoring a con¬ 
stant term that does not affect the dynamics, this gives 

Hab = HojaA^A + HcobB^B, (3) 

with normal mode eigenfrequencies 

I ± a/(A2 - - i6G^A(jjyn 

^A,B = \ -^-. (4) 


They are plotted in Fig. 1 as a function of the detuning A. It 
is straightforward to see that for A «: -oJm, the polariton “A” 
is photon-like and the polariton “B” is phonon-like, while in 
the opposite limit -ojm A < 0, it is the polariton “A” that is 
phonon-like and the polariton “B” is photon-like. 

Consider then a situation where the phonon reservoir is at 
some finite temperature Tphonon, while the optical field is cou¬ 
pled to a reservoir at T = 0 - an excellent approximation 
at visible frequencies - and concentrate first on the polariton 
“B” only. It is possible to realize an Otto cycle for that nor¬ 
mal mode in the following way : Start from the system in 
thermal equilibrium at a large negative detuning A, in which 
case “B” is essentially at the temperature of the phonon bath, 
with corresponding thermal excitation number = Nb, 

and adiabatically change A across the resonance A = -oJm 
and toward small negative value close to 0 (to avoid the onset 
of instabilities). In that first adiabatic stroke the polariton “B” 
changes its character from phonon-like to photon-like, and the 
energy of the thermal phonons is converted into intracavity 
photons that perform work on the oscillating mirror via radi¬ 
ation pressure - more on that in the following section. Once 
the detuning has reached its final value the system is then al¬ 
lowed to thermalize with the cavity field reservoir at T = 0 
(first thermalization stroke), releasing heat in the process. The 
third stroke is again adiabatic (second adiabatic stroke). It 
consists in changing A back to its large negative value. Fi¬ 
nally the cycle is closed by allowing the polariton “B”, which 
has now regained its phonon-like character, to thermalize with 
its reservoir at Tphonon by absorbing heat (second thermaliza¬ 
tion stroke.) These four strokes are sketched schematically in 
Fig. 1(b). 

The polariton ‘A” simultaneously also goes through a ther¬ 
modynamic cycle, with however significant differences. First, 
it is initially coupled to a reservoir at T 0, so that the ini¬ 
tial thermal polariton occupation is Na ~ 0. Second, the first 
thermalization stroke for the “B” cycle, which takes a time 
of the order of a few is not long enough to also thermal¬ 
ize the ‘A” polariton provided that the optical damping rate is 
much faster than the mechanical damping rate, /c » y, which 
is normally the case in optomechanical systems. Under these 
conditions the population of mode ‘A” remains essentially un¬ 
changed and equal to zero, and the ‘A” cycle does not produce 
any work (positive or negative.) That is, provided that the 
changes in detuning A can be realized in a perfectly (quantum) 
adiabatic fashion, the two cycles remain completely decou¬ 
pled. This, however, no longer holds if quantum adiabaticity 
cannot be maintained. This will have important consequences 
in the context of the dispersive quantum measurements of sec¬ 
tions IV and V. 


III. OUTPUT WORK 

We now turn to a discussion of the work performed by the 
optomechanical heat engine. It is useful to carry out this dis¬ 
cussion both in the normal mode picture and the bare modes 
picture, if only because the measurements we have in mind 
will be on the optical field. 
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(a) 




FIG. 1: (Color online) (a) Frequencies of the two normal modes (po- 
laritons) of the optomechanical system for G/du^ = 0.1 in the red- 
detuned case A < 0. The dashed lines correspond to the of the bare 
photon and phonon modes. The plot also indicates that for large neg¬ 
ative detunings the “B” polariton is phonon-like, and photon like for 
small negative detunings, (b) Sketch of the four strokes of the Otto 
cycle for the “B” polariton. See text for details. 


For notational convenience we decompose the energy oper¬ 
ator formally as 

U = W+Q. (5) 

Although there is still considerable debate in the literature 
about the nature of a work operator, as mentioned in the In¬ 
troduction, this issue needs not concern us here as we are only 
concerned about operational ways to evaluate its first moment 
and fiuctuations in a specific experimental setting. 

The first law of thermodynamics, cast in infinitesimal form 
can then be formally written as 

d{U) = d{W) + d{Q), (6) 

where and {W}, are the mean energy, heat, and work, 

respectively. For a general open quantum system with density 
operator p and Hamiltonian H we have 

{U) = Tvipm, (7) 

so that in the Schrodinger representation the average (clas- 
sical) values of the infinitesimal work and heat increments 
are lld^ 


d{W) = Tr[p{dH)\, (8) 

d{Q) = Tr[{dp)H]. (9) 

Our main focus is the measurement of the work produced 
by the system. For an isolated quantum system, that work 
may be determined unambiguously via a two-time measure¬ 
ment process However this approach is problematic 

for an open quantum system since it would require additional 
measurements on the reservoirs il. We therefore adopt an 
operational approach based on stochastic quantum thermody¬ 
namics m. 

Specifically, our starting point is the description of the evo¬ 
lution of open quantum systems in terms of a large ensem¬ 
ble of N quantum trajectories that are solutions of a 

stochastic Schrodinger equation of the general form 

d\i//(t)) = (Ddt 1^dw)\i//(t)), (10) 


where the superoperator Ddt accounts for both the Hamil¬ 
tonian evolution of the system and non-unitary contributions 
that account for dissipation and decoherence mechanisms as¬ 
sociated with measurement processes, and the stochastic term 
Hdw, where dw describes one or more Wiener process of 
zero mean with dw^ = dt, accounts for the stochastic quan¬ 
tum jumps resulting from both reservoir noise and/or quantum 
measurement The initial conditions \i//j(0)) are ran¬ 

domly selected consistently with the initial density operator 
of the system. 

The properties of any quantum observable can be calculated 
based upon the statistics deriving from these trajectories. In 
particular for each trajectory, we can compute the work by 
integrating the Eq. ([5]) as 

X tf 

{llfj{t)\ — \i^j{t))dt ( 11 ) 

where ti and tf are the initial and final times. The resulting 
mean value and the variance of the work are then obtained as 


^ TI/ 

Z W j 

7=1 


( 12 ) 


and 


, ^ (Wj - {W})^ ,, , , 

^ ^ J ' = {W^) - {wf. (13) 

7=1 


N 


In the limit N ^ oo, the statistics resulting from the quantum 
trajectories approach the correct result. The numerical simu¬ 
lations presented in section V are based on that approach. 


A. Output work in the optomechanical QHE 

We have seen that in the normal mode picture the “B” 
polariton heat engine is driven through an Otto cycle. Pro¬ 
vided that dissipation is weak enough to be negligible dur¬ 
ing the adiabatic strokes only work is performed during those 
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strokes, while heat is only exchanged during the thermaliza- 
tion steps El. Since the “B” polariton population is Nb = 0 
after thermalization with the optical heat bath at T = 0 we can 
then restrict the determination of the output work to the first 
adiabatic stroke, where 

d{Q) = Tr[(dpAB)HAB] = 0, (14) 

and 

d{W} = TrlpABidHAB)] = NbMojb. (15) 

Here we have used the polariton Hamiltonian ([3]) and assumed 
that Na ~ 0, as previously discussed. Since in the adiabatic 
stoke Nb is conserved, the average work is simply given by 
the change in energy of mode “B”, 

{iy) = ^ ' dw = (16) 

where A/ and A/ are the initial and final detunings, see Eq. ([T]). 
Note also that in the case of perfect adiabaticity, the “B” po¬ 
lariton number distribution gives directly the full statistical 
distribution of the work performed by the system as well. 

In the bare mode representation, the first stroke is not adi¬ 
abatic, and the photon and phonon distributions are time de¬ 
pendent. Equation ([8]) then takes the form 

d{W) = Tvlpatmab)] = (17) 

where we used the Hamiltonian (O and ha - {a^d) is the av¬ 
erage number of excitations in the photon mode a. In this 
picture the work is then given by 


length 

A(y) = Ao - (19) 

where y is a classically controlled length change, assumed 
small with respect to the total cavity length, gM = g/ym is 
the optomechanical coupling normalized to the mirror zero- 
point motion , and Aq is the initial detuning. With Eq. (O 
we can then express the work in terms of the spatial integral 
of the position-dependent radiation pressure force, 

r F^{y)dy, (20) 

dji 

where 

^rp(y) = Hunaiy)- (21) 

A possible classical scheme to measure the work output of 
the engine is illlustrated in Eig. [2l The optomechanical res¬ 
onator comprises the oscillating end mirror driven by radi¬ 
ation pressure and a mirror of large mass M whose classical 
position can be controlled externally by the potential y(y) pro¬ 
vided by a piezoelectric element, therby varying the detuning 
A(y) in the presence of the radiation force Fr,. To use a ther¬ 
modynamical metaphor, we may think of the input mirror as 
a classical piston that is pushed by the expanding photon gas. 

The total system Hamiltonian is then 

HT = Hab + HM, ( 22 ) 


{W} = -n j ha(A)dA. (18) 

If the stroke is perfectly adiabatic the values of the average 
work obtained from expressions (O and (fTSl) are equal. How¬ 
ever if either the optical or the mechanical damping is signifi¬ 
cant on the time scale of the adiabatic stroke, or if the variation 
of the optical detuning induces non-adiabatic transitions and 
in particular a non-vanishing population of polariton A, then 
Eq. (O is no longer exact. The expression of the average 
work (fTSl) in terms of the mean photon number remains how¬ 
ever valid. 


where Hat is given by Eq. (O and 


(23) 


is the classical Hamiltonian for the massive control mirror. 
The classical equations of motion for that mirror are then 


dt 

dp 

dt 


OHt 
OHt 


P_ 

M’ 


dv(y) 

■ ^y 




(24) 

(25) 


IV. QHE WORK MEASUREMENT 

We now discuss several possible measurement schemes that 
can be considered to quantify the work performed by the heat 
engine and its fiuctuations. To set the stage we first consider 
a simple classical approach before considering two types of 
quantum measurements. 


where Hr is the classical limit of the total Hamiltonian Hr. If 
M is large enough that it can be considered as infinite com¬ 
pared to all other optomechanical elements we have dy/dt ^ 
dpidt ^ 0. That is, the force exerted by the control system 
balances the expectation value of the radiation pressure force. 


dV(y) _ 
dy " ■T” 


(26) 


A. Classical measurement scheme 

We can think of implementing the variation in detuning A 
required for the adiabatic stroke through a change in cavity 


This shows that provided the kinetic energy of the large mirror 
remains essentially zero all work performed by the photons is 
converted to the control potential energy and can be measured 
in that way. 






5 



FIG. 2: (Color online) Schematic setup for the classical measurement 
of the output work, whereby the work performed by the radiation 
pressure acting on the mirror of large mass M can be stored in the 
control system. See text for details. 


B. Continuous quantum measurements 

It is not possible to directly monitor the occupation of the 
polariton mode “B” since it consists of quasiparticles that are 
coherent superpositions of photon and phonon states. Instead, 
and in analogy with the classical measurement scheme we ap¬ 
ply a weak continuous measurement scheme El to monitor 
the intracavity photon number and calculate the total work by 
performing the integral in Eq. (fTSl) . Similarly to the classical 
scheme the total system Hamiltonian is then of the form 

Ht = Hat + Vm, (27) 

with Vfn describing the weak interaction between the photon 
and the measuring quantum system. Since in general does 
not commute with the optomechanical Hamiltonian Hat the 
measurements lead in general to back-action on the QHE that 
affects output work and its efficiency. In the following we 
will use the method of quantum trajectories to simulate the 
measurement processes and investigate their infiuence on the 
mean work (W) and its fiuctuation AW^. 

Operationally we consider continuous measurements of the 
intracavity field realized by passing through the resonator a 
dilute beam of two-level atoms that interact weakly with it, 
with at most one atom at a time inside the resonator. The 
state of the field is then inferred from a measurement on the 
atoms after they exit the cavity. We study both the cases of 
absorptive and dispersive atom-field inteactions, see Eig.[3l 


1. Absorptive measurements 



FIG. 3: (Color online) Schematic setup for a continuous quantum 
measurement of the output work with a beam of two-level atoms, a) 
Absorptive measurement: the cavity field is resonant with the atomic 
transition and the coupling induces real oscillations in the atomic 
population, which result in the loss of intracavity photons, b) Dis¬ 
persive measurement: the cavity mode frequency is far off-resonant 
from the two-level atom transition frequency, resulting in a disper¬ 
sive interaction that only modifies the phase of the atomic ground 
state wave function. See text for details. 


surements on the reduced density operator for the optical field 
is then governed by the stochastic master equation \4^ 

dp(t) = -Aa {2dpd^ - d^dp - pd^d^ dt 

^[Aa [dp + pd- 2{d + d^)p^ dw (29) 
where dw is a Wiener process and 

4 = sIt (30) 


is a measure of the strength of the measurement, r being the 
transit time of an individual atom through the resonator. Note 
that in obtaining Eq. ([29l) the time increment dt is assumed 
to be long compared to the atomic transit time r, so that this 
equation describes the statistical effect on the field of a large 
number of atomic measurements. 

The first term on the right-hand side of Eq. (l29l) accounts 
for the additional dissipation channel of the intracavity field 
resulting from absorption by the successive atoms, and the 
second term describes the stochastic changes on the intracav¬ 
ity field about its expected value (a -i- as a result of the 
measurement outcomes. 


Consider first the resonant situation where the atomic tran¬ 
sition frequency ojeg = oJc and the atoms are prepared in their 
ground state \g} before being injected inside the optical cavity, 
see Eig. 3(a). The atom-field coupling is given in the rotating 
wave approximation by 

Vm = nga(a^&ge + O^eg) (28) 

where ga is the single-photon Rabi frequency of the transition 
and d-fj = j = {e, g}. The effect of the continuous mea¬ 


2. Dispersive measurements 

We now turn to the situation where the interaction between 
the two-level atoms and the intracavity field mode is off- 
resonant. Upon adiabatic elimination of the upper electronic 
stateit is described by the effective Hamiltonian 

Vm = figda^ai&ee - <Tgg) = figda^a{&+- + &-+), (31) 
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which conserves the mean photon number fia = {a^a). Here 
gd = the off-resonant coupling between the intensity 

of the field and the energy of the atomic levels. 

The atoms are now prepared in the superposition |-r) = 
(k) + k))/ V 2 and information on the intracavity field is in¬ 
ferred from a change in phase of the atomic state. In that 
situation the effect of the measurements on the optical field is 
described by the stochastic master equation 

dp{t) = ^Ad(2hapha - hip - phl)dt 

+ yfddinaP + pha - 2{ha)p)dw, (32) 

where Ad = glr. 

As was the case for resonant coupling this equation also 
comprises two contributions, the second one accounting for 
the stochastic changes of the mean intracavity intensity (h) 
about its expected value resulting from successive measure¬ 
ments. But because of the quantum non-demolition nature of 
the non-resonant atom-field interaction, the dissipative chan¬ 
nel of Eq. (I 29 I) is now replaced by a number conserving term 
that results in an additional damping of the phase of the optical 
field. 

Importantly, in the specific case of our optomechanical 
QHE the effective interaction (ISTt couples the two polariton 
branches “A” and “B”, transferring excitations between these 
two modes. As we see in section V this can be thought of as 
a nonadiabatic coupling that has in general a significant effect 
on the work that can be extracted in the Otto cycle of the “B” 
polariton . This is in contrast to absorptive measurements, 
as the interaction (1^ . does not significantly couple the two 
normal mode branches. Still, both measurement schemes re¬ 
sult in the appearance of additional photon loss channels that 
limit the amount of extractable work. In both cases these mea¬ 
surement back-action mechanisms may be viewed as heat ex¬ 
change between the engine and the environment. 


V. NUMERICAL RESULTS 


This section presents selected results from numerical simu¬ 
lations of the continuous measurement of the work output of 
the QHE and its fiuctuations as defined in equations ([T^ and 
(US, both for absorptive and dispersive measurements. The 
numerical results were obtained by averaging for each choice 
of parameters 20,000 trajectories obtained from a combina¬ 
tion of the Hamiltonian evolution of the system as the detun¬ 
ing A(0 is varied across the Otto cycle and the solution of the 
stochastic Schrodinger equations 


dm = 


1 2 

(fifl - (na)) 


dt 


and 


dm = 


+ {na})dw]\i//(t)), (33) 


--Aaid^d - {d + d^)d + ^ ' ) 


dt 


+ ^[Ta{d - |(A(0), (34) 

corresponding to the stochastic master equations ([29l) and (13^ 
to account for the continuous measurements. 

To guarantee that the adiabatic strokes are indeed adiabatic 
in the absence of measurements, it is important to change A(0 
sufficiently slowly that non-adiabatic transitions between the 
two polariton branches remain negligible, but fast enough that 
the damping of both the optical field and the mechanical os¬ 
cillator at rates k and y respectively, remain negligible. This 
is particularly the case near the avoided crossing at A = -cOm- 
The insert of Eig. [4] shows as an example the time evolution 
of A(0 (in units of cOm) used in the simulations of the first 
adiabatic stroke to minimize this problem, with A(0 changing 
rapidly away from the avoided crossing and very slowly in its 
vicinity. As a result non-adiabatic transitions remain negli¬ 
gible, as illustrated in the black curves (labeled by a square) 
of Eig. [H The population of the “B” polariton mode remains 
essentially constant during the adiabatic stroke, and the ‘A” 
polariton population remains essentially zero. 

The additional curves in Eig. [4] show the time dependence 
of the average populations ^a( 0 and Nsit) during the first 
stroke of the heat engine, again neglecting mechanical and 
optical damping, for both dispersive (red lines with triangles) 
and absorptive (blue lines with circles) measurements. (Note 
that Na remains extremely small during that stroke and its evo¬ 
lution is nearly indistinguishable from the situation without 
measurements.) This is the most important stroke as far as 
extracting work from the engine is concerned, since the pop¬ 
ulation of normal mode “B” remains extremely small during 
the second adiabatic stroke as a result of its thermalization at 
the optical reservoir temperature T = 0. In these examples the 
temperature Tphonon of the phonon bath and the initial detuning 
A(0) are chosen such that Nb{0) ^ 4 and Na(0) ^ 0. 

Comparing the evolution of the mean populations of the po¬ 
lariton modes for absorptive and dispersive measurements il¬ 
lustrates clearly the difference in their back-action on the op¬ 
eration of the QHE. In the case of absorptive measurements 
the “B” polariton population decreases significantly during 
what would otherwise be an adiabatic, population-conserving 
stroke. Since that stroke occurs fast compared to that 
damping results solely from an additional photon dissipation 
at rate Aa that results from the measurements, see Eq. (l29l) . 
Importantly, though, these measurements do not result in any 
significant transfer of population to the ‘A” polariton. 

The situation is qualitatively quite different for dispersive 
measurements. In that case there is no significant loss in total 
polariton population, but instead a significant transfer of popu¬ 
lation from mode “B” to mode “A”. This is because the effect 
of dispersive measurements is an additional source of deco¬ 
herence, but no loss of population, see Eq. ([3^ . Dispersive 
measurements change the frequency of the photons stochas¬ 
tically, as seen by the term proportional to d'^d in the master 
equation ([^ . and thereby they change the structure of the po- 
laritons. In contrast, in the absorptive case the measurements 
remove excitations from the system, but without affecting the 
structure of the polaritons. 

The population transfer between the normal modes “B” 
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FIG. 4: (Color online) Time evolution of the mean excitations Nb and 
Na of the polariton modes “B” (solid lines) and “A” (dashed lines) 
during the first stroke of the heat engine, averaged over 20,000 trajec¬ 
tories of the stochastic Schrodinger equation. Black lines marked by 
squares: no measurement. Red lines marked by triangles: Dispersive 
measurement with Ad = O.Oddu^. Blue lines with circles: absorptive 
measurements with Aa = O.Oddu^- Other parameters: G = 

A/ = A/ = -OAojm ,/c = 5 x and y = I0~^cx)fn. Inset: 

time dependence of the pump-cavity detuning A(t) in units of 
Time in units of l/du^- 


and “A” associated with dispersive measurements causes a 
reduction in the work performed by the system that can be¬ 
come quite dramatic due to the resulting unavoidable cou¬ 
pling between the two normal modes. Since for the polariton 
branch ‘A” the Otto cycle is reversed and produces positive 
work EH), that is, work is performed by the environment on 
the polariton one can even reach situations where the 
effective available work of the two systems, which are inex¬ 
tricably coupled, becomes positive. For absorptive measure¬ 
ment, in contrast, the two normal modes remain essentially 
uncoupled, and although the output work can be significantly 
reduced it always remains negative. 

Figure Oa) shows the mean value of the output work -{W) 
for increasing measurement strength Aa,d, illustrating its re¬ 
duction due to measurement back-action. Surprisingly per¬ 
haps for equal measurement strengths dispersive measure¬ 
ments cause a stronger reduction in work than absorptive mea¬ 
surements. Since the thermalization processes are not affected 
by the measurement scheme the heat absorbed by the system 
from the mechanical reservoir, (Qm), remains the same for all 
scenarios. For this reason, the efficiency of the quantum heat 
engine follows directly from the work as 


(- 1 ^) 

(Qin) 


(35) 


where the minus sign accounts for the fact that by convention 
the work done by the engine is negative ((W) < 0). 

We now turn to the fluctuations of the output work. Their 
variance, plotted in Fig. [3b), shows a signiflcant quantitative 
difference between the situations for absorptive and disper¬ 
sive measurements. In the first case (dotted line) the fluctu¬ 
ations decrease monotonically as a function of the measure¬ 



FIG. 5: (Color online) (a) Expectation value and (b) variance of the 
output work , in units of for the full cycle as a function of 
the measurement strenght (in units of ojm). The points are results of 
numerical simulations and the lines serve to guid the eye. In both fig¬ 
ures, the sqaures (red solid line) and circles (blue dashed line) stand 
for the dispersive and absorptive measurement scheme respectively. 
The statistic is done upon 20000 trajectories. The stroke times are 
ti = 40du;;^ t = ^3 = 406 c >>;;^ ^4 = 4 x All other 

parameters as in Fig.|4| 


ment strength, while in the case of dispersive measurements 
they remain roughly constant. 

One can gain a better understanding of this behavior from 
the probability distribution P(-W) of the output work as a 
function of measurement strength. Figure [6] shows this dis¬ 
tributions in the absence of measurements and for two mea¬ 
surement strengths, for both dispersive and absorptive mea¬ 
surements. Without measurements the probability distribution 
consists of a series of discrete peaks that correspond to the dis¬ 
tribution of Fock states in the initial thermal distribution of the 
mechanical oscillator. Assuming perfect adiabaticity each of 
these mechanical Fock states is converted into a photonic Fock 
state during the first adiabatic stroke, and produces a specific 
amount of work. (The width of the peaks is due to residual 
non-adiabatic effects.) 

Continuous measurements result in a broadening of the 
peaks, an effect of the stochastic nature of the detection pro¬ 
cess (see upper panels of Fig. [6]) and, in the case of absorptive 
measurements, a decrease in amplitude of all peaks except the 
one corresponding to the vacuum field, a consequence of the 
additional photon decay channel. This is the reason for the re¬ 
duction in variance as the measurement strength is increased. 
In contrast, for dispersive measurements the distribution shifts 
toward positive values of the work. This is more apparent 
in the lower panels, which show the same distribution on a 
logarithmic scale. This is a direct consequence of the cou¬ 
pling with the ‘A”-polariton engine which, as we have seen, 
tends to be characterized by positive work. Because absorp- 
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live measurements don’t couple the polariton modes in any 
significant way, this effect is almost completely absent in that 
case. Finally, since the mean total number of polaritons in 
modes “A” and “B” varies only slightly over the chosen mea¬ 
surement strengths, and the changes in the photon distribution 
are much less significant than for absorptive measurements, 
resulting is weak changes in the variance of the extracted work 
as a function of measurement strength. 

VI. SUMMARY AND OUTLOOK 

Summarizing, we have developped a measurement model 
to characterize the mean work and its fluctuations in an op¬ 
tomechanical QHE and performed a numerical study of the ef¬ 
fect of continuous quantum measurements on its performance. 
We considered measurement schemes involving the continu¬ 
ous monitoring of the intracavity photon field, with both dis¬ 
persive and an absorptive interactions with a dilute beam of 
two-level atoms. By determining the average value and the 
variance of the work we are able to quantify the measurement 
back-action effects. In both cases, the measurements were 
found to induce a reduction in the average work performed 
by the engine and thus a reduction in its efficiency. However, 
the detailed reasons behind these reductions are qualitatively 
different. In the dispersive regime, the measurement induces 
transitions between the two polariton modes, and hence two 
thermodynamic cycles, one producing negative and the other 
positive work. The final result is a reduction of the efficiency 


with increased fiuctuations in the work output. In the absorp¬ 
tive detection scenario, in contrast, photons are lost from the 
system via the interaction with the quantum probe that acts as 
an effective (zero temperature) reservoir. In this case, both the 
average value and the fiuctuations of the work decay mono- 
tonically. 

Perhaps the most intriguing result of this study is the real¬ 
ization that quantum measurements permit to control the op¬ 
eration of coupled QHE, and to switch the operation from a 
cycle transferring energy from a hot to a cold reservoir to the 
reverse situation. This suggests that it might be possible to 
consider quantum heat pumps whose operation is controlled 
by the back-action of quantum measurements. This and other 
aspects of QHE will be further explored in future work. 
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